Finite-time Singularities in Surface-Diffusion Instabilities are Cured by Plasticity 



OO 
O 

o 
in 

(N 



I 

C/3 



I 

c 

o 
o 



> 
^- 
o 

m 

(N 
O 
OO 

O 



Ting-Slick Lo. Anna Pomyalov, Itamar Procaccia and Jacques Zylbcrg 
-Depi 0/ Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

A free material surface which supports surface diffusion becomes unstable when put under ex- 
ternal non-hydrostatic stress. Since the chemical potential on a stressed surface is larger inside an 
indentation, small shape fluctuations develop because material preferentially diffuses out of inden- 
tations. When the bulk of the material is purely elastic one expects this instability to run into a 
finite-time cusp singularity. It is shown here that this singularity is cured by plastic effects in the 
material, turning the singular solution to a regular crack. 



We address the problem of an amorphous material with 
a free surface on which the material can diffuse such that 
the surface normal velocity is proportional to d 2 p,/db 2 
where \x is the local chemical potential and b is the pa- 
rameterization of the interface. When the material it- 
self is purely elastic, this phenomenon leads to an insta- 
bility which was termed 'thermal grooving' by Mullins 
[l[ who discovered it. The phenomenon of grooving is 
equally present in crystals, where it provides an annoy- 
ing mechanism for the failure of growing crystals 0, El , 
as it does in a range of amorphous solids that concern us 
here. Mullins considered the linear instability taking into 
account only the curvature dependence of the chemical 
potential. In fact, the chemical potential along the inter- 
face is strongly dependent on the elastic energy; linear 
stability analysis taking both effects into account [H re- 
veals that the surface is stable for short wavelengths but 
unstable for longer ones, with a usual 'fastest growing 
mode' whose wavelength depends on the material pa- 
rameters. The late stage of development of this insta- 
bility was initially studied by a handful of researchers, 
namely Asaro and Tiller [B[ and Grinfeld Q and followed 
by many 0, SHU- The result is that the instabil- 
ity runs into a finite-time singularity, with the growing 
indentation forming a cusp. Clearly, this often explored 
0, floL 11 1 mathematical phenomenon cannot be physical, 
and its discovery leaves open the question of the physical 
mechanism that may cure the singularity. 

The question of what might cure the finite-time sin- 
gularity in the Asaro-Tiller- Grinfeld (ATG) instability 
remained dormant until recently Brener and Spatschck 
proposed that inertial effect in the velocity of the moving 
boundary may tame the singularity Q. These authors 
pointed out that without inertial effects the velocity of 
the tip v appears in one dimcnsionless combintation, i.e. 
VTq/D where r$ is the radius of the tip and D the diffu- 
sion coefficient (of dimension length 4 /time, and cf. Eq. 
(jT5j)). Therefore there is no mechanism to select v or 
ro, and as rrj decreases without limit, v increases with- 
out imit. Once inertial effects are taken into account the 
velocity appears also in the combination v/vr where vr 
is the Rayleigh wave speed. Thus a selection of both v 
and tq can hapen. While clearly correct, the present au- 
thors stress that in many cases the surface diffusion is 



very slow, leading to small interface velocities which do 
not justify the incorporation of inertial terms. We fo- 
cus here on such cases where the question of taming the 
cusp-singularities remains open. 

In this Letter we propose that the generic mechanism 
for the taming of the ATG instability may be plastic de- 
formation in the stressed material, especially near the pu- 
tative cusp. To test and demonstrate this proposition we 
will employ the recently proposed theory of elasto-plastic 
dynamics in amorphous systems 12|. To this theory, 



which is valid in the bulk of the material, we couple the 
surface diffusion, allowing the chemical potential to take 
its stress dependence from the elasto-plastic theory. For 
concreteness we choose to explore this interesting physics 
on the inner surface of a circular hole which is stressed at 
infinity in a radial fashion. The surface diffusion modifies 
the shape of the slightly perturbed circular hole, leading 
eventually to a highly non-linear morphology. With a 
sharpening interface due to the surface diffusion insta- 
bility, stresses in the bulk increase rapidly, exceeding at 
some point in time the yield stress of the material, trig- 
gering plastic flows which are dissipated by the exertion 
of plastic work 12|, 13, (3 ■ ^ is interesting to observe the 



coupling of both processes, namely surface diffusion and 
plasticity, as they become competitive and of opposite 
influence on the morphology, to a point where the finite 
time singularity is removed. In addition to shedding new 
light on the late stage of the ATG instability we find that 
the elasto-plastic theory employed here, which is sensi- 
tive to the plastic properties of matter, allows a natural 
understanding of this a-priori seemingly hard problem. 

The model system that we consider here is an infinite 
2-dimcnsional isotropic elasto-plastic sheet with a hole in 
the center whose radius is R(0). For r(9) < R(9) the sys- 
tem is void, whereas the elasto-plastic material occupies 
the region r(6) > R{6) . The boundary is traction free, 
meaning that on the boundary cryHj — where er is the 
stress tensor and n is the unit normal vector. The equa- 
tions of acceleration and continuity are exact, reading: 
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Here the full material derivative T> is defined for an arbi- 
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trary tensor A as: where u> is the spin tensor = § f |j£ - Reading 

Eq. ([T]) in radial coordinates we get: 
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Here s and r arc defined via the transformations 
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Below our velocities arc sufficiently small to allow ne- 
glecting the nonlinear terms —v 2 /r and v r vg/r. On the 
other hand nonlinear terms containing derivatives are re- 
tained, since the derivatives arc large. 

The velocity at the interface, R(0), reads: 



— = -dgR + v r 



(7) 



When the surface evolves, the stresses in the bulk 
evolve accordingly. A fundamental assumption of our 
elasto-plastic theory is that the total rate of deformation 
D tot = i[Vi) + (Vd)'] can be a represented as a linear 
combination of its elastic and plastic components 



D> 



D e 



D pI 



(8) 



Here the elastic contribution, D el , is assumed to be lin- 
early dependent on the stress (linear elasticity) 
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where K and [i are the 2-dimcnsional bulk and shear 
moduli and p and Sij are the pressure and the deviatoric 
stress tensor, respectively. The plastic rate of deforma- 
tion, D pl , is determined by a set of internal fields which 
are discussed at length in 12| where the elasto-plastic 



theory is presented in detail. For the purpose of this Let- 
ter it is enough to state that the tensorial field m acts 
as a 'back-stress' due to plastic deformations, and the 
scalar field x 1S the effective temperature that controls the 
amount of configurational disorder in the elasto-plastic 
materials. The constitutive relations that were derived 
for these fields read 
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Vx/Vt = e xT(s ij ,mi j )(xo -x) 
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In these equations all the stresses were normalized by the 

yield stress of the material s y , using s = SijSij /2Sy. 

The function C(s) has been chosen to make the material 
relatively brittle. 

To this theory we need to couple now the surface dif- 
fusion, expressed in terms of the normal velocity, v n (0) 
on the boundary. Without the effects of elasto-plasticity 
in the bulk the normal velocity satisfies 



Z^ 2 ^ d 2 n 
k B T db 2 



(15) 



where D s , the surface diffusion constant, 0, particle vol- 
ume and 5, the number of particles per unit area. In 
solving the coupled problem the total normal velocity 
should be computed as a sum of this contribution and 
the one coming from Eqs. (U]) and |5]). 

The chemical potential on the boundary, fi(0), is asso- 
ciated on the one hand with the destabilizing curvature 
and on the other hand with the stabilizing surface energy 



H = Ho — -fK + £ 



(16) 



Here is the chemical potential of the unperturbed sur- 
face, the curvature k(0) = (R 2 + 2R' 2 - RR'^/CR 2 + 
R' 2 ) 3 / 2 and the strain energy density £ = ^crijaj [l5j. 7 
is the surface energy. 

In terms of these effects on the chemical potential one 
derives the equation for the normal velocity due to sur- 
face diffusion alone: 
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where v is the Poisson ratio and E Young's Modulus. 
After non-dimensionalization and projecting from v n to 
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v r one ends up with the following equation in terms of 
dimcnsionlcss quantities (denoted with tilde): 
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As noted, we need to couple Eqs. (JH) , © and (fT8]) to be 
solved together, such that the surface normal velocity is 
made from the sum of contributions coming from the bulk 
dynamics and the surface dynamics respectively. This 
should be done while keeping the traction free boundary 
conditions and the initial condition of a pure elastic so- 
lution of a slightly perturbed circle. In practice we used 
the fact that the elastic response is the fastest process 
in this problem. Accordingly we solved at each itera- 
tion first the elastic part of the model (with D pl = 0) to 
find the stress fields which is in agreement with the given 
interface, without taking into account any plastic defor- 
mation. Second, elasto-plastic relaxation was allowed to 
take place, until the system reached elasto-plastic equi- 
librium, allowing the interface to change. Here 'equilib- 
rium' means that D pl is smaller than 10 -4 . Lastly, a step 
of surface diffusion was allowed to take place using an 
adaptive time step such as to bound the maximal move- 
ment of the boundary by 10~ 4 . The last step changes 
the morphology of the boundary again, necessitating a 
re-calculation of the elastic fields around the new bound- 
ary, etc. Since the effect on the velocity of the interface in 
the last two steps is additive, these steps (being infintes- 
imal) could be also done simultaneously with impunity. 
We chose to separate the last two steps since the plastic 
and surface diffusion processes arc non-dimcnsionalized 
independantly and have different normalization values of 
characteristic times and stresses. 

In order to realize an infinite sheet it is convenient to 
transform the (r, 9) coordinate system through a confor- 
mal transformation to a new, finite domain (£,#) with 
Ce [0,1]: 



C(9)=R(9,t)/r 



(19) 



In the finite space all the derivatives are computed using 
finite differences and are redefined for the transformed 
space using the chain rule: d Xi = d Xi Xkd Xk . Explicitly 

(20) 

At £ = all the derivatives vanish and on the boundary 
C = 1 the time derivatives are estimated by linear ex- 
trapolation in the £ direction. For the sake of numerical 
stability small viscosity terms were added to the accel- 
eration components. The coupled equations were solved 
using K = 100 and /i = 50. Since the speed of sound is 
orders of magnitude larger than the velocity of the inter- 
face, we could safely neglect the effects of Eq. ©, giving 




FIG. 1: A typical half profile of the stressed interface under 
the action of surface diffusion and plastiticity. To the bare 
eye the effect of plasticity is not seen here, and one needs to 
compare elastic and plastic solution in Figs. [5] and [3] 
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FIG. 2: The tip curvatures of the elastic solution. This solu- 
tion appears to approach a singularity at t* = 5.37 x 10 -4 . 



up on seeing the sound waves. The initial condition on 
the perturbed circle were 



R(t = 0,9) = 1 + 0.01 cos 29 



(21) 



the stress at infinity was chosen to be = 0.9s y . The 
Poisson ratio v in Eq. (fl~8"|) is 1/3. The surface diffu- 
sion equations contain fourth order derivatives, calling 
for spectral techniques for sufficiently stable evaluation. 
All the other derivatives were computed by finite differ- 
ences. 

Results and Discussion: The typical morphology of 
the unstable interface is shown in Fig. [1] The elastic 
solution for the finite-time singularity is well established 
and was reproduced in our numerics. The curvature k 
in the growing cusp first grows exponentially and rapidly 
switches to faster regime that agrees with the growth law 
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To make this growth obvious we plotted the curvature of 
the elastic solution in Fig. [2] as a function of (t* — i)~ 1//2 
with t* = 5.37 x 10~ 4 . Once plasticity is allowed to inter- 
vene, it prevents the finite-time singularity by blunting 
the tip and by dissipating the stress. 
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FIG. 3: The tip velocity (panel A), its first and second time 
derivatives (panels B and C respectively). We see that the 
singularity is cured and the velocity decelerates due to the 
plastic effects. 



On the boundary, the smoothening of the interface in 
the vicinity of the cusp via blunting is the "cure" of the 
singularity. The ever-increasing curvature occurring in 
the elastic solution is prevented in the plastic solution 
by the plastic flow induced by stress concentration. The 
avoidance of the singularity is shown by the deceleration 
in the tip velocity, see Fig. [3] panels B and C. 

It is important to stress that although the plasticity in 
the bulk succeeds to cure the finite-time cusp-singularity, 
the role of surface diffusion is far from being negligible. 



Without it, the stressed circle would remain stable to 
small shape fluctuations, as was demonstrated recently in 
[l6| . The surface diffusion makes the circle unstable, and 
the instability results in the growth of a groove. With- 
out plasticity in the bulk the solution loses its meaning 
at t = t*, whereas now, with plasticity playing its use- 
ful role, the solutions continue to exist at times t > t*, 
in a form of a lengthening groove, or crack, whose tip is 
protected from cusping by the plastic effects. At some 
point the crack will increase its velocity due to the Grif- 
fith mechanism, and then the problem becomes inertial 
again. 
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